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We study pairing of an interacting three component Fermi gas in two dimensions. By using a 
mean field theory to decouple the interactions between different pairs of Fermi components, we study 
the free energy landscapes as a function of various system parameters including chemical potentials, 
binding energies, and temperature. We find that the s-wave pairing channel is determined by both 
chemical potentials and the interaction strengths between the three available channels. We find a 
second order thermal phase transition and a series of first order quantum phase transitions for a 
homogenous system as we change the parameters. In particular, for symmetric parameters, we find 
the simultaneous existence of three superfluid orders as well as re-entrant quantum phase transitions 
as we tune the parameters. 



I. INTRODUCTION 

Recent experimental progress achieved in ultra-cold 
atomic gases allows one to set up test beds for con- 
trolled study of many body physics. By tuning three 
dimensional two-body scattering length between atoms 
in different hypcrfine spin states of a dilute system at low 
temperatures, interaction strength between atoms can be 
controlled very precisely [1J. This can be done by using 
magnetically tuned Feshbach resonance |2J. Moreover, 
interaction and spatial dimensionality can be effectively 
controlled by applying an optical lattice. An effective 
two-dimensional system can be created by applying a rel- 
atively strong one-dimensional optical lattice to an ordi- 
nary three-dimensional system. As there are always two 
body bound states exist for attractive two body poten- 
tials in two dimensions [3] , different pairs of Fermi atoms 
can undergo Bose Einstein condensation and form super- 
fluidity at low temperatures. The two dimensional bound 
state energies can be controlled by tuning either three di- 
mensional scattering length or the laser intensity which 
used to create one dimensional lattice to accommodate 
2D layers. 

In this paper we study three component Fermi gases in 
two dimensions. A mixture of 6 Li atoms which has favor- 
able collisional properties among its lowest three hyper- 
fine spin states will be an ideal system to explore novel 
three component superfluidity. Three component 6 Li 
mixture in three dimensions have already been trapped 
and manipulated experimentally [H |S] . In Ref . |3] , us- 
ing radio frequency spectroscopic data and a quantum 
scattering model, scattering lengths and the Feshbach 
resonance positions in the lowest three channels of 6 Li 
atoms have been determined. As there are three broad 
s-wave Feshbach resonances, one can prepare the system 
at various interaction strengths between each pairs. In 
Ref. [5], collisional stability of the lowest three channels 
of 6 Li atoms has been studied. As the spin relaxation 
time is large compared to the other time scales in the 
experiments, experimentalists were able to maintain a 
fixed spin population throughout the experiments. As 



a physically accessible system, a three-component ultra- 
cold atomic system can be used to study the physics of 
nuclear matter. Three-component Fermi pairing is be- 
lieved to occur in the interior of neutron stars and in 
heavy-ion collisions ;6]. Nevertheless, this system can be 
used to understand the competition between quantum 
phases and re-entrant phase transitions. 

Properties of three-component Fermi systems have 
been extensively studied in recent past [Tj. I5| f9|. [TO ] fTT I fTS] . 
However, all these studies are carried out in a three di- 
mensional or one dimensional environments. Further, 
authors in all these references except Ref. [5] have re- 
stricted their parameter space either by assuming equal 
interaction strengths or equal chemical potentials or by 
neglecting the interaction between some hyperfine spin 
components. In ref. [H], the authors have studied the 
properties of a harmonically trapped three component 
gas in three dimensions. 

In this paper, we neglect the possibility of three body 
bound states in two dimensions and consider only two 
body pairing states. This is reasonable, because of the 
system we are considering is dilute and the atomic in- 
teractions are short range in nature. Therefore, it is un- 
likely to have many atoms interacting in the same region 
of space. Further, we neglect the harmonic confinement 
and consider the system as spatially homogenous. Two- 
component Fermi gases in two dimensions have already 
been studied in theory jT3j [14] . The purpose of this paper 
is to investigate how pairing will occur when a third spin 
component is added to such a two-component gas. More 
precisely, we study the competition of individual com- 
ponents to form Bose condensed pairs by investigating 
the landscapes of the free energy as a function of vari- 
ous parameters which include the temperature, chemical 
potentials of the hyperfine spin components, and the in- 
teraction strengths between different pairs of fcrmions 
components. For appropriately chosen parameters, we 
find that the system undergoes a second order thermal 
phase transition from normal state to a superfluid state 
as one lowers the temperature. At low temperatures, 
we find a series of first order quantum phase transitions 
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FIG. 1: Free energy contours showing a second order thermal 
phase transition for the parameters fj.2 = Q-8(ii, [1$ = 0.75/ii, 
Ebi = 0.1/ii, i?s2 = jUi and _Eg3 = 0.99/ii. From (a) to (d) 
temperature varies as ksT — /ii/1.00, fcsT = fii/l.lljksT = 
/ii/1.20, and fcsT = /^i/50. The global minimas in figures (a) 
to (d) are at (Ai/^i, A 2 /mi, A 3 /mi) = (0,0,0), (0, 0, 0.23), 
(0, 0, 0.79), and (0, 0, 1.66) respectively. 



where r 2 — x 2 + y 2 , V2D is the 2D gradient operator and 
U nn i is the 2D interaction strength between component 
n and n' . The operator ipn{ r ) creates a fermion of mass 
m with hyperfine spin n = 1,2, 3 at position r = (x,y). 
The chemical potential of the n'th component is \i n . No- 
tice that we have neglected the interaction between the 
same components. This is reasonable as we are consid- 
ering a dilute atomic system, and the interactions are 
short-range in nature, s-wave scattering channel is dom- 
inated over the other scattering channels. By using a 
mean field decoupling of the interacting terms, the mean 
field Hamiltonian in the momentum space can be written 
as, 

Hmf = £ tiAirfs + \ S>Ji?«^ + h.c) - £ ^(2) 

where i — k,n and j — —k,n. Here we defined the 
superfluid order parameters Uij(ipiipj) — Ay and two 
matrices A and B, 



as we change the chemical potentials or the interactions 
between hyperfine spin components. At low tempera- 
tures, we find simultaneous existence of three types of 
superfluid phases (at symmetric parameters correspond- 
ing to different pairing channels) and normal phases in 
this novel three-component Fermi system. Interestingly, 
we find that the system can undergo re-entrant phase 
transitions as we simultaneously tune the chemical po- 
tentials and interactions. 

The paper is organized as follows. In the following sec- 
tion, we introduce the theoretical model and use a mean 
field approximation to decouple the interaction terms. 
Then using a canonical transformation, we diagonalize 
the Hamiltonian to derive the free energy of the system. 
In section III, we present our results with a discussion. 
Finally, our summary and conclusions are given in section 
IV. 



II. FORMALISM 

We consider an interacting three-component Fermi 
atomic gas trapped in two dimensions. We take the 
model Hamiltonian of the system as 



(ex \ 
A = e 2 (3) 
V e 3 / 



/ A 3 -A 2 \ 
B = I -A3 Ax (4) 
V A 2 -Ai / 

where e n — h 2 k 2 /(2m) — [i n and Ay = eij k A k . As the 
mean field Hamiltonian is quadratic in Fermi operators, 
it can be diagonalized with a canonical transformation to 
get 

H MF = £ AnvUn + \ £(e„ - K) ~ £ ^ (5) 

n : k n,k n 

where we use the notation J7y = \eij k \U k . The Fermi op- 
erators r/ n represent the quasi particles in the system with 
n = 1,2,3. The quasi particle energies A„ = y/Xn are 
given by A„ = 2^Q cos{[ 6> + (n - 1)2tt]/3} - A k /3. The 
parameter 9 = arccos[i?/ \J— Q 3 ] with Q = (3B k — A 2 )/9 
and R = (9A k B k - 27C k ~ 2A 3 k )/54 [15]. Here we defined 
Ak = - En fn - 2 £ A* , B k = £ n A^ + 2 £„ A 2 n e 2 n + 

En^m A„A m + 1/2 J2n^m e n e m + J2n^m^l ^n e mQ and 



FIG. 2: Superfluid order parameters A3 as a function of 
temperature. We fixed the binding energies Ebi = 0.1//1, 
Eb2 = 0.5/ti, Eb3 = fJ-i and chemical potentials /i2 = 0.8/ti 
and ^3 = 0.75/ii. 
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where /3 = l/(k B T) is the inverse temperature and fcs 
is the Boltzmann constant. As the short range nature of 
the interaction, the grand potential is diverging so that 
regularization must be done in standards way by writing 
-\A n \ 2 /U n = J2 k \A n \ 2 /(h 2 k 2 /m + E Bn ). Here E Bij = 
\iijk\EBk is the binding energy between two hyperfine 
spin components i and j. Notice that we use the same 
notation for Es n as we used for A^- = e^^Afc and Uij — 
\tijk\Uk- For two dimensions, converting the J^ fc into 
integral J d 2 k/(2-K) 2 and then by changing the variable 
by k 2 — z, the grand potential can be converted into an 
one dimensional integral. We numerically perform this 
integral and numerically minimize the grand potential 
for the seven parameter space (three chemical potentials, 
three binding energies and the temperature). 



III. RESULTS AND DISCUSSION 

In Fig. [TJ we plot the free energy landscapes for dif- 
ferent temperatures at a selected set of parameters. We 
choose the binding energy in the 2 — 3 channel to be 
small (Ebi = 0.1/ii) so that pairing is not possible in 
this channel. As a result, we find a second order thermal 
phase transition as we lower the temperature. As can 
be seen, at high temperature [Fig. [Q(a)], free energy is 
minimum when both paring order parameters (A2 and 
A3) in channels 1—3 and 1—2 are zero. As we lower the 




FIG. 3: Low temperature free energy contours showing first 
order quantum phase transitions for the parameters /ii = 
M2 = = Mi Ebi = O.l/i, Eb3 = 0.6/i and fcsT = ju/50. 
From panel (a) to (c) Eb2 varies as O.6/1, 0.7/j, re- 

spectively. In panel (a) and (c) the minimas are given 
at (Ai/fi, Az/fi, A 3 ///) = (0,0,1.25) and (0,1.37,0) respec- 
tively. In panel (b), free energy gives many stationary points 
where the global minimum is a quarter of a circle. In panel 
(d), we use Ebi = Eb3 = 0.05/i and the minimum is given at 
(0.46, 0, 0). 



temperature, channel 1 — 2 undergoes pairing and form 
Bose condensation. This is because the binding energy 
and average chemical potential in this channel is larger 
than those of channel 1 — 3 and 2 — 3. Further lower- 
ing the temperature results more atom pairing and con- 
densation in channel 1 — 2 giving larger superfluid order 
parameter A3. In principle, it is possible to have a se- 
quence of second order thermal and first order quantum 
phase transitions at three different critical temperatures, 
if one change the interactions or the chemical potentials 
together with temperature. The reason for this sequence 
of phase transition is that there are many ways of pairing 
when various favorable channels are available. 

In Fig. [2j we plot the temperature dependence of the 
superfluid order parameter (A3) in channel 1 — 2 for cho- 
sen values of parameters. We choose the parameters such 
that pairing is possible only in channel 1 — 2 so that a 
single minimum is available in the free energy. As can 
be seen, the superfluid order parameter continuously in- 
creases as one lower the temperature, showing a second 
order thermal phase transition. 

By varying the average chemical potentials and bind- 
ing energies of the pairing channels at low temperatures, 
one can control the first order quantum phase transitions 



from one superfluid phase to another. As a demonstra- 
tion, we plot the free energy landscapes in Fig. [3] for a 
selected set of parameters. Again, we chose the parame- 
ters such that the pairing in channel 2 — 3 is very weak 
and the corresponding superfluid order parameter (Ai) is 
zero. When the binding energy in channel 1 — 2 is larger 
than that of the channel 1 — 3, the free energy minimum is 
at a non zero value of A3 but zero value of A2 at the same 
chemical potentials. However, when the binding energies 
are equal at equal average chemical potentials, both su- 
perfluid order parameters are non zero and free energy 
gives many stable stationary points as seen in Fig. (3[(b) 
(now the global minimum is not a point, but a quarter of 
a circle). The reason for this line of global minimum is 
the symmetry of the parameter space. By increasing the 
binding energy in channel 1 — 3 over the channel 1 — 2, the 
minimum free energy pass to the non zero A2 but zero 
A3. Similar first order quantum phase transitions can 
be seen by controlling the average chemical potentials at 
fixed and equal binding energies. More generally, by con- 
trolling the chemical potentials and binding energies, one 
can observe not only a series of quantum phase transi- 
tion, but a phase with multi-component superfluid order 
(simultaneous existence of three superfluid order param- 
eters). Notice that one can have a re-entrant quantum 
phase transition by changing both chemical potentials 
and binding energies simultaneously. 

In Fig. [4] we plot superfluid order parameters as a 
function of the chemical potential of the third compo- 
nent. The binding energies are fixed to be the same for 
all three channels. As can be seen in figure, fcrmions 
pairing occurs in channel 1 — 2 at smaller 113. This is 
because the average chemical potential in this channel is 
the largest for /13 < Eb- Further, as the average chemi- 
cal potential is constant, superfluid order parameter A 3 
is constant. For /13 > Eb, average chemical potential in 
channel 2 — 3 is the largest and increasing with increasing 
/i3. As a result, superfluid order parameter Ai increases 
with fj,3. In the entire range of /i3, average chemical po- 
tential in channel 1 — 3 is smaller than that of the other 
channels so that the pairing in this channel is not favor- 
able. As seen in Fig. |4j the superfluid order parameters 
have large discontinuity which represents a sharp first 
order quantum phase transition. 

At the same chemical potentials and the same in- 
teraction strengths of the channels, free energy gives 
many stable stationary points at which the condition 
Af + A| + A§ = C is satisfied. The constant C de- 
pends on both the chemical potentials and the interac- 
tion strengths (binding energies). As shown in Fig. [5] 
when the free energy has a minimum, the order parame- 
ters represent a surface in order parameter space. 

Within our mean field description, we were able to 
handle only two-body correlations. One needs to go be- 
yond mean field theory to understand the role of three- 
body correlations in a three-component system. If three 
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FIG. 4: Superfluid order parameters Ai/Eb (black line) and 
A^/Eb (gray line) as a function of [izjEs- We fixed the bind- 
ing energies Ebi = Ebi = Ebz = Eb and the temperature 
fcgT = Eb /50. The chemical potentials are /ii = Eb and 
fi2 — \.2Eb- In the entire range of (13, A2 is zero. 



atoms can overlap in the same region of space, then the 
three-body correlations can play a role giving Thomas ef- 
fect [T7] and Efimov effect [H]. Thomas effect is the col- 
lapse of a three body system due to the overlap of atoms. 
Atom loss in a trap is undoubtedly related to the Thomas 
effect. In a three component atomic system, collision be- 
tween a condensed pair and a third species atom can sup- 
port the Thomas effect. Efimov effect is the accumulation 
of three-body bound states at strongly interacting limit. 
These effects are forbidden in two component gases. For 
sufficiently low densities, these effects are forbidden even 
in three component systems so that our results are appli- 
cable to dilute ultra-cold atomic gases. We discussed the 
pure 2D limit in this paper, however one can generalize 
the theory to include the weak atom tunneling between 
layers as done in Ref. [H] for two component gases. 

Our results in two dimensions look qualitatively simi- 
lar to the ones obtained in three dimensions in Ref. [7]- 
However, we find that the superfluidity is more sensitive 
to the parameters in two dimensions than three dimen- 
sional systems. As we have seen above, superfluidity is 
very sensitive to both chemical potentials and the inter- 
actions between different pairs. In current experimen- 
tal setups, Feshbach resonance allows one to control the 
interactions to different values by tuning the scattering 
lengths. However, the scattering lengths between differ- 
ent pairs of fermions cannot be controlled independently 
Therefore, chemical potential is the suitable parameter to 
drive the quantum phase transitions. Typically, chemi- 
cal potentials can be controlled by changing the atomic 
population in both two dimensions and three dimensions. 
However, in order to change the chemical potentials this 
way, one has to start the experiment all over with a dif- 
ferent atomic sample. The advantage of using quasi two 
dimensional system is that one can change the effective 
chemical potential by controlling the tunneling between 
layers. This tunneling can be controlled by the laser in- 
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FIG. 5: Superfluid order parameters for symmetric param- 
eters. We use the same chemical potentials for the three 
species /ii = /i2 = M3 = M an( A same interaction strengths 
for the three channels Ebi = Eb2 = Eb3 = Eb- The three 
surfaces (A? + A| + A§ = constant) shown in the order pa- 
rameter space are for Eb = 0.5/^, 1.2/Lt, and 2.0/^. We fix the 
temperature to be ksT — p/50. 

tensity of the optical lattice. As shown in Ref. [TS], the 
first order quantum phase transitions in two dimensional 
systems can easily be controlled by the laser intensity. 



IV. SUMMARY AND CONCLUSIONS 

We studied Fermi superfluidity of an interacting three 
component system in two dimensions. We used a mean 
field theory to investigate the behavior of free energy as a 
function of chemical potentials, binding energies, and the 
temperature. Depending on the chemical potentials and 
binding energies, we find a second order thermal phase 
transition as we lower the temperature. At low tempera- 
tures, by controlling the parameters, on can have a series 
of first order quantum phase transitions and re-entrant 
phase transitions. These series of phase transitions are 
associated with pairing between different hyperfine spin 
components of fermions. 

The possible pairing is determined by both average 
chemicals potentials and the interaction strengths of the 
paring channels. The channel which has largest paring 
strength forms the superfluid, while the unpaired com- 
ponent form a Fermi sea. At low temperature, first order 



quantum phase transition can be induced by increasing 
the average chemical potential or the interaction strength 
of one channel over the other. If the pairing strengths 
are equal in all three channels, then it is possible to have 
three superfluid phases simultaneously. As we have not 
considered the interaction between condensed pairs, we 
do not expect the phase separation of superfluid phases 
in spatially homogenous environments [9]. However in 
trapped systems, it has been shown that the phase sep- 
aration of superfluid phases is possible in three dimen- 
sions |8J. 

We speculate that the simultaneous existence of three 
types of superfluid phases in trapped systems can be de- 
tected by standard experimental methods. For exam- 
ple, superfluidity can be demonstrated by the creation of 
vortices [T5] and then distinguished them by probes cou- 
pling to each atom types. Alternatively, one can measure 
the energy gap using radio frequency spectroscopy [5D] 
or measure condensate fraction using the pair projection 
method [21]. 
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